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ABSTRACT 

Aims. We study the grain properties and location of the forsterite crystals in the circumstellar environment of the pre-planetary nebula 
(PPN) IRAS 17150-3224 in order to learn more about the as yet poorly understood evolutionary phase prior to the PPN. 

Methods. We use the best-fit model for IRAS 17150-3224 of |Meixner et al.| ( [2002^ and add forsterite to this model. We investigate 
different spatial distributions and grain sizes of the forsterite crystals in the circumstellar environment. We compare the spectral bands 
of forsterite in the mid-infrared and at 69 pm in radiative transport models to those in ISO-SWS and Herschel/PACS observations. 
Results. We can reproduce the non-detection of the mid-infrared bands and the detection of the 69 pm feature with models where the 
forsterite is distributed in the whole outflow, in the superwind region, or in the AGB-wind region emitted previous to the superwind, but 
we cannot discriminate between these three models. To reproduce the observed spectral bands with these three models, the forsterite 
crystals need to be dominated by a grain size population of 2 pm up to 6 pm. We also tested models where the forsterite is located in 
a torus region or where it is concentrated in the equatorial plane, in a disk-like fashion. These models show either absorption features 
that are too strong or a 69 ^m band that is too weak, respectively, so we exclude these cases. We observe a blue shoulder on the 69 
pm band that cannot be explained by forsterite and we suggest a possible population of micron-sized ortho-enstatite grains. 

We hypothesise that the large forsterite crystals were formed after the superwind phase of IRAS 17150-3224, where the star developed 
an as yet unknown hyperwind with an extremely high mass-loss rate (>10“^ Mo/yr). The high densities of such a hyperwind could 
be responsible for the efficient grain growth of both amorphous and crystalline dust in the outflow. Several mechanisms are discussed 
that might explain the lower-limit of ~2 ^m found for the forsterite grains, but none are satisfactory. Among the mechanisms explored 
is a possible selection effect due to radiation pressure based on photon scattering on micron-sized grains. 

Key words. Stars: evolution. Stars: AGB and post-AGB, Infrared: stars. Stars: mass-loss, Techniques: spectroscopic, Stars: individ¬ 
ual: IRAS 17150-3224 


1. Introduction 

IRAS 17150-3224 (also known as AFGL 6815 or 
OH 353.84-1-2.98) is a pre-planetary nebula (PPN), charac¬ 
terised by a stellar remnant surrounded by a previously expelled 
dust and gas envelope. The infrared emission from the dust en¬ 
velope of IRAS 17150-3224 has been studied by |Guertler et al.| 
( |1996| l, who modelled the spectral energy distribution (SED) 
of that source up to wavelengths of 1.3 mm. Using a spherical 
geometry for the dust envelope they found that they needed 
either large grains or a second cold shell in order to properly 
fit the far-infrared emission of IRAS 17150-3224. Based on 
Hubble resolved images, |Ueta et ffiT] ( |2000| l and 
( 1999| l later showed that the envelope of IRAS 
axisymmetric. |Meixner et ar] ( |2002| l modelled both the SED and 
the resolved images using an axisymmetric dust distribution and 
amorphous grains with a minimum grain size of 0.001 pm and a 
maximum size of up to sizes of several hundred micrometres. 


Meixner et al.| 
7150-3224 is 


The axisymmetric envelope of IRAS 17150-3224 is the ma¬ 
terial lost by IRAS 17150-3224 during its asymptotic giant 
branch (AGB) phase. The axisymmetric morphologies are also 
recognised in planetary nebulae themselves and these morpholo¬ 
gies are thought to predate the PPN phase. 


IRAS 17150-3224 is likely an evolutionary product of a 
massive (>5 Mg) main-sequence star ( [Meixner et al.|20()2) . This 
links IRAS 17150-3224 to stars that in their AGB phase stay 
oxygen-rich. Such massive AGB stars have mass-l oss rates of the 
order of ~10 Mg/yr when in their OH/IR phase (Vassiliadis & 
Wood|[T99^. This period of high mass-loss is often called the 


superwind (jlben & Renzini|1981|[Knapp & Morris|1985[[Bedijn[ 
[TM7l[Woodetal.|1992|l. 


Even though the superwind of massive AGB stars was intro¬ 
duced so that the central star could lose enough mass to reach 
white dwarf masses, it is now becoming increasingly clear that 
the superwind is not capable of shedding enough mass. The su- 
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perwind timescale turns out to be too short (<2000 yrs) and the 
number of superwind phases seems to be limited to only one. 
These short timescales have been shown with several methods, 
either based on gas and SEP modelling ([Heske et al.|1990[|Just-| 


tanont et al. 2006||2013[[Decin et al.|2007HJusttanont et al.|i996[ 

Chesneau et al.|2005[ Groenewegen|2012[) or based on the anal- 

ysis of forsterite spectral features ( de Vries et al.|20f4] l. 

These short timescales pose an evolutionary problem. 
The mass loss rates of massive AGB stars can go up to 
15 ■ 10“^ Mo/yr and even slightly higher (Justtanont et al.|l996 


Justtanont & Tielens|1992[|Schutte & Tielens|1989HGroenewe-^ 

gen 1994| l. Such mass-loss rates, combined with the found 
timescales, gives a total mass lost in one superwind of ~Q.2- 
0.6 Mq (depending on the gas over dust ratio of 100-300, |Lom- 
[baert et al.|2013[ ). Massive AGB stars need to lose several solar 
masses before they can leave the AGB. Therefore one would ex¬ 
pect to see several superwind phases as extended shells around 
massive AGB stars and their remnants, in the same way as ex¬ 
tended shells are seen around carbon-rich AGB stars ( |Cox et ah] 
2012 Maercker et al.|20f^. It is intriguing that even though sev¬ 


eral oxygen-rich AGB stars have been observed with Herschel, 
none of them show any extended structure ( Cox et al.|2012| l. It 
is likely that this is an indication that massive AGB stars have 
no extended shells and thus no previously ejected superwinds. 
It is then hard to understand how these massive AGB stars lose 
enough mass to evolve away from the AGB into a white dwarf 
de Vries et aI.|(20T4ll hypothesise an intermittent phase be¬ 


tween the PPN phase and the superwind, where the star has a 
mass-loss rate one or two orders higher than during the super¬ 
wind. This phase of extreme mass-loss could be referred to as 
the hyperwind. With hyperwind mass-loss rates of the order of 
10“^-10“^ Mo/yr, several hundreds to thousands of years would 
be enough for the central star to evolve into a white dwarf. 

The central star of a massive AGB during its superwind is 
already completely obscured and its SEP is very red, peaking 
at ~30 pm. The SEP of a star in its hyperwind phase would be 
even redder and would peak at even longer wavelengths, while 
its central star would still be totally obscured. So far, stars in 
their hyperwind phase have not been observed or identified, be¬ 
cause the phase is so short, the superwind is still too optically 
thick to see the hyperwind through it, or such sources might 
have been wrongly classified (possibly as protostellar). How¬ 
ever hints of the hyperwind can already be seen in the mod¬ 
els of IRAS 17150-3224 and other post-AGB sources (e.g. 
I RAS 16432-3814 in|Pijkstra et al.|2003]or IR AS 18276-1431 
in |Murakawa et al.|2013| l. |Meixner et al.| ( 2002 [ ) show that, at its 


highest, the mass-loss rate of IRAS 17150-3224 must have been 
of the order of 8.5-10-3 Mo/yr. This makes IRAS 17150-3224 
a crucial object to study in the light of this evolutionary conun¬ 
drum. 

Another property that links IRAS 17150-3224 to massive 
AGB stars is the presence of forsterite (Mg 2 Si 04 ) in its outflow. 
Massive AGB stars in their superwind phase show prominent 


forsterite features at among others 11.3, 33.6 and 69 pm (Waters 

et aL|1996 

[Molster et al.|20021[de Vries et al. [20101 [Blommaert 

eraL120M 

1 [de Vries et al.[[2014[l. The abundance of forsterite 


in these superwinds can be as high as 14% by mass (de Vries 


|erar][20T0l [20T^ . Since the spectrum of IRAS 17150-3224 
does not have any convincing mid-infrared spectral features of 
crystalline dust species (Meixner et al.|(2002|), see Eig.[T]), it 


thought that the outflow of IRAS 17150-3224 contains no crys¬ 
talline dust component, but recently Blommaert et al.|((2014 1 re¬ 


ported the detection of the far-infrared spectral band at 69 pm of 
forsterite. 
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Fig. 1. The main window shows the ISO-SWS spectrum of 
IRAS 17150-3224 in grey l |Sloan et al.||2003y The continuum con¬ 
structed over the 9.7 pm band of amorphous olivine (as explained in 
sect. |4.3^ is shown as a dashed black line. The top-left inset shows the 
optical depth of the 9.7 pm absorption feature (see sect. |4.3^ in grey and 
in dashed black the continuum under the 11.3 pm band of forsterite. 
In the bottom right inset we show the 69 pm band as observed with 
Herschel/PACS in grey ( [Blommaert et al.]|2014| l and the constructed 
continuum under the 69 pm band in dashed black. 


The study of the sharp spectral bands of forsterite is of great 
value since these spectral bands are dependent on many proper¬ 
ties of the grains (see sect. [3 1 and thus on the outflow in general. 
As we will show in this work, a special combination of parame¬ 
ters is needed for IRAS 17150-3224 in order to explain both the 
absence of forsterite features in the mid-infrared and a strong de¬ 
tection at 69 pm. To this end we use radiative transport models. 

In this paper we start by introducing the observations in 
sect. |2] The behaviour of the different mid- and far-infrared 
bands of forsterite are described in sect. [3 Sect. |4] contains the 
modelling strategy of this work, the radiative-transport code we 
use and how we study the spectral features in the model and ob¬ 
served spectra. We continue in sect.|^with the results and end the 
paper with a discussion in sect.|^and our conclusions in sect.j^ 


2. Observations 


Most recently pVleixner et al.| ( |2002] l have modelled the SED and 
images of IRAS 17150-3224. In their fitting they have included 


the B-band Hubble Space Telescope image from Ueta et al. 
( |2000| l, the 9.8 pm images from |Meixner et al/] ( | 1999| ) and the op¬ 
tical and infrared photometry of a compilation of data made by 
Ueta et al. ( 2000|l (including data from van der Veen et al. ( |1989] l. 


Hu et al.| ( |1993| l, [Kwok et^ ( |1996| ), Reddy & Parathasarathy 


(1996), Meixner et a l. (1999), Ueta et al. (2000), and IRAS). 
Meixner et al.| ( |2002[ l were able to reproduce the SED and ir 


ages of IRAS 17150-3224 with an axisymmetric density dis¬ 
tribution and by including a population of amorphous olivine 
grains with a size distribution from 0.001 pm up to ~200 pm 
and then falling off exponentially. 

We will build on the results of [Meixner et al.| ( |2002| l by 
studying the now known presence of forsterite in the outflow of 
IRAS 17150-3224. Eor this we use the ISO-SWS spectrum of 
IRAS 17150-3224 ( [Sloan et al.|2003 see Eig.[T]|. The ISO-SWS 


spectrum shows a very smooth dust continuum and a 9.7 pm ab¬ 
sorption band, which are both due to amorphous silicates. As 
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mentioned by |Meixner et al.| ( |2002| l the ISO-SWS spectrum con¬ 
tains no clear indications of any crystalline silicate bands (for 
example at 11.3 jum or 33.6 /rm). The presence of three weak 
emission features at 33.6 jum, ~38 fim, and ~41 jum could be ar¬ 
gued, but this is uncertain at the end of the wavelength range of 
ISO-SWS. The fact that an absorption band of amorphous sili¬ 
cates is seen at 9.7 yum while no absorption band of forsterite is 
seen at 11.3 /tm is strikingly different from that seen in obser¬ 
vations of massive AGB stars, the likely predecessor phase of 
IRAS 17150-3224. 

In contrast to the absence of mid-infrared features, the far- 
infrared spectrum of IRAS 17150-3224 shows a clear 69 jum 
band of forsterite ( [Blommaert et al.|2014| see Fig.[T]l. In the far- 
infrared we use the Herschel/PACS observations of the 69 jum 
(Groenewegen et al. 201 1[ Blommaert et al.|2014|) band. The de- 


tection of the 69 jim band in IRAS 17150-3224 was reported 
by Blommaert et al. ( 2014| l to be broader than could be ex¬ 
plained from its central wavelength position. [Blommaert et al. 


( |2014| l showed that a temperature gradient for the forsterite com¬ 
ponent was not enough to explain the width of the 69 yum band 
of IRAS 17150-3224. In sect. we will explain in more detail 
how the band properties of the 69 fim band depend on grain tem¬ 
perature and grain size. 

In this work we will study the combination of the non¬ 
detection of features in the mid-infrared together with the de¬ 
tection of the 69 yum band. We will also show that the broadness 
of the 69 yum band as reported by jBlommaert et ar] ( |2014| l is due 
to a shoulder on the blue side of the 69 yum band, which cannot 
be due to forsterite. 


3. Forsterite features 

As many different investigations have shown, the spectral fea- 


2003 1 [Suto et al. [2006| [Bowey et al. 

|2002| [Sturm et al. 

2010 

Mulders et al. 

20111 

de Vries et al. 

20121 Sturm et al. 

2013 

Maaskant et a 

1. 2015 

[Blommaert et al. 20141 de Vries et al. 


whose shape is strongly dependent on the grain temperature, the 
iron content of the crystal and, for micron-sized grains, also on 
the crystal size. The band broadens and its central wavelength 
position shifts to the red as either the temperature increases, 
the iron content increases or the grain size becomes larger than 
7 jum in size ([Koike et al. 2003 [ Suto et al.|2006t [Stum et al.[ 


2013[ Maaskant et al.|2015 i. Fig. ^ shows the behaviour of the 

69 yum band as a function of grain size. It can be seen that the 
effect of grain size on the 69 fj.m band is small below grain sizes 
of 7 yum, while above this limit the 69 fim band broadens signif¬ 
icantly as a function of grain size. 

The bands of forsterite in the region of the 9.7 yum band of 
amorphous silicates, like the 11.3 yum band, originate from de¬ 
grees of freedom within the Si 04 tetrahedral and they are there¬ 
for not sensitive to the composition of the crystalline olivine 
( Koike et al.|2003 1 . By approximation, these bands are also not 
sensitive to the grain temperature ( [Zeidler|2012[ l. The forsterite 
bands in the mid-infrared (like for example the 33.6 yum band, 
see Fig. are, compared to the 69 yum band, mildly dependent 
on the grain temperature and composition ([Koike et al.[[2003[ 
[Zeidler|2012| l. 

All spectral bands of forsterite are dependent on grain size, 
but at which grain size their dependence becomes significant de¬ 
pends on the wavelength position of the band. The grain size 
dependence of the bands comes from the fact that at a certain 
size, the interior of the grain is not sampled by the radiation field 



10 15 20 25 30 35 



68.0 68.5 69.0 69.5 70.0 70.5 71.0 

Wavelength (^nll) 


Fig. 2. The effect of grain size on the mid-infrared and 69 yum bands of 
150 K forsterite grains. The grey, blue, green, red, teal, purple and yel¬ 
low curves are for forsterite with grain sizes of 0.001, 0.1, 1, 3, 5, 7 and 
9 pm, respectively. The opacities are calculated for grains with a DHS 
shape distribution ( [Min et al.|2005[ l, using the optical constants of [Suto[ 
[et al.| ( |^06| l. The dashed black line are the opacities when calculated 
with CDE (Continuous Distribution of Ellipsoids, [Bohren & Huffman[ 
[!^ 


anymore, causing the band to broaden and eventually to disap¬ 
pear. As is shown in Fig. already at grain sizes of 1 pm the 
band at 11.3 pm becomes weaker as a function of grain size. At 
a grain size of 1 pm this is not yet the case, causing the band at 
that wavelength to still be quite insensitive to grain size. Even¬ 
tually the 69 pm band will be influenced by the grain size when 
the grain size becomes larger than ~1 pm. 

Because of the grain size dependence, the peak ratios of 
the spectral bands are a good probe of the grain size. This 
was recently shown and used by [Maaskant et alT] ( [2015[ l. They 
found ~10 pm forsterite crystals in the proto-planetary disk of 
HD 141569 by studying the relative strengths of the 69 pm band 
compared to mid-infrared bands. However, before doing such a 
grain size analysis, a good model for the temperature and density 
structure of the circumstellar environment needs to be known. 
For IRAS 17150-3224, a satisfactory model was derived by 


Meixner et al. (20021, which provides us with the needed pa¬ 


rameters to calculate the optical depth, density and temperature 
of the envelope in order to study the forsterite spectral signature 
(see sect. E^- 


4. Method 

4.1. Strategy 

Here we discuss three ways by which the spectrum of forsterite 
dust would only show a 69 pm band and no features in the mid- 
infrared. The first is that the material is very cold and the tem¬ 
perature of the material has a black-body spectrum that peaks at 
the far-infrared wavelengths or slightly beyond. A second option 
is that the medium has an optical depth that causes all the mid- 
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infrared features to be exactly in between absorption and emis¬ 
sion. The third option is the presence of micron-sized forsterite 
crystals (>1 pm), which do not show spectral bands at short 
wavelengths because the grain size is significantly larger than 
the wavelength, but still have the 69 pm band in the far-infrared. 

In order to find out which of the three options is the cause 
for the peculiar forsterite spectrum of IRAS 17150-3224, we 
study forsterite features over a wide range of wavelengths, from 
11.3 pm up to the 69 pm band. We choose features situated at 
wavelengths far apart, which helps us to probe the largest possi¬ 
ble ranges in optical depth and temperature. The 11.3 pm band 
is very sensitive to the grain size in the range between 0.1 pm 
and 5.0 pm, but not to the grain temperature and crystal compo¬ 
sition (iron content of the crystalline olivine grains). The 69 pm 
band is not sensitive to the grain size (at grain sizes <7 pm), but 


is very sensitive to the grain temperature (for example see Suto 


et al. 120061 IKoike et al.|2003l|Sturm et al.|20101|de Vries et al. 

2012 


The density structure of the dust around IRAS 17150-3224 
was determined by |Meixner et al.| (|2002|l (see sect. |4.2| and Ta¬ 
ble We take their derived parameters for our radiative trans¬ 
port models and add forsterite to the model (see Table |^. The 
presence of forsterite in the outflow has no effect on the density 
or temperature structure and the results of |Meixner et al.| ( [2002| l, 
since the opacities of forsterite do not add any continuum opac¬ 
ity to the dust and only a small amount is needed to reproduce 
the 69 pm band (see Table |^. 

In order to extensively test all possible cases that could ex¬ 
plain the features of IRAS 17150-3224, we explore five distribu¬ 
tions for the forsterite, schematically depicted in the left column 
of Fig. 0 


- Model Whole: forsterite follows the same density distribu¬ 
tion as the amorphous dust and is present in both the super¬ 
wind and the low mass-loss rate AGB-wind previous to the 
superwind. 

- Model Superwind: forsterite is only present in the superwind 
(R < Rsw)- 

- Model AGB-wind: forsterite is only present in the outflow 
that happened previous to the superwind (R > Rsw)- 

- Model Torus: forsterite is only present in the densest part of 
the superwind (R < Rsw. G > 50°) 

- Model Disk: forsterite is only present in a disk-like region in 
the superwind (R < Rsw, 0 > 85°) 


Here 0 is the angle from the pole and R is the radial distance 
from the centre of the star. For all models, in the region where 
forsterite is present in the outflow, it follows the density distribu¬ 
tion of the amorphous silicates (see Eq.[^. 

For these five models we test a grid of six grain sizes for 
forsterite: <0.1, 1.0, 2.0, 4.0, 6.0 and 8.0 pm. We do not consider 
smaller grain sizes for forsterite since the absorption opacities do 
not change when the size is decreased from 0.1 pm (see Fig.|^. 
For these different grain sizes we run a grid of models with dif¬ 
ferent forsterite abundances from 0.5% up to 25% (for the Disk 
model only did we increase the abundance to 90% in order to get 
a signal at 69 pm). In these grids we search for the models with 
abundances that correctly predict the strength of the observed 
69 pm band. The abundances found for the different models are 
listed in Table Among the models that reproduce the 69 pm 
band we search for one that also explains the absence of the mid- 
infrared bands in the spectrum of IRAS 17150-3224. 
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Table 1. Best-fit model parameters of Meixner et al. i 2002 1 


Meixner et al. 

20021) fit parameters 

Lstar ( Lq) 

27,200 

Tstar (K) 

5200 

Rstar (R©) 

201 

Distance (kpc) 

3600 

Rin (AU) 

648.4 

Rsw (AU) 

1944.4 

Rout (AU) 

64 ■ 10^ 

Vexp (km/sec) 

15 

Inclination 

82.0° 


Total dust mass (Mq) 

-^meix 

Bmeix 

c 

meix 

Emeix 

F 

^meix 

F 

^ meix 


Am. grain size 
Am. grain size slope 


4.6 ■ 10-2 
159 
2.0 
1.5 
1.0 
4.0 
1.5 

0.001 - 200 pm 
-3.5 


Table 2. Forsterite abundances used for the models in Fig.[^ 


Forsterite abundance (%, see sect.|4.l| 


Grain size (pm): 

<0.1 

1 

2 

4 

6 

8 

Whole 

1.3 

1.3 

1.3 

1.6 

1.8 

2.1 

Superwind 

3.8 

3.8 

3.8 

3.8 

3.8 

3.8 

AGB-wind 

1.6 

1.6 

1.6 

2.6 

3.4 

4.8 

Torus 

10.6 

10.6 

10.6 

11.5 

12.3 

11.9 

Disk 

90.0 

- 

- 

- 

- 

- 


4.2. Models 


[Meixner et al.|p002] l fitted the SED up to 1.1 mm as well as the 
optical and near-IR resolved images of IRAS 17150-3224. It is 
beyond the scope of this work to improve on this already satis¬ 
factory model, and we adopt its parameters for the central star 
and the circumstellar environment (see Table [1]). In Sect. |4.4| 
we do test the robustness of our results when uncertainties are 
introduced in the total dust mass and the stellar effective temper¬ 
ature determined by [Meixner et al.| ( 2002| l. [Meixner et al.| ( |2002| l 
concluded that the best fit was obtained with an equatorially en¬ 
hanced density distribution. 


p{R, 0) = Pin I ^ 

\K\n 


-fi(l+Csin''(6»)(e-^^^^sw) /e-^^in/Rsw) )) 


( 1 ) 


x[l +A(l -cos(0))''(e-<^/’2sw)‘^/e-''^'”^'2™)°)]. 


with Rin < R < Rsw; R is the distance from the centre of the 
central star and 0 is the angle from the polar direction (0 - 0); 
Rin and pin are the inner radius and the density at the inner radius 
of the dust shell; and Rsw is the superwind radius outside of 
which the outflow is just spherically symmetric. [Meixner et al.[ 
( |2002| l create a time-independent outflow (in all directions) by 
using an r^^ density distribution (B = 2.0). 


The parameters of the best-fit of [Meixner et al.[ ( |2002[l are 
shown in Ta ble [T] The total dust mass determined by Meixner 


et al. 


(2002 I is 4.6 ■ 10 Mq. This mass is constrained by the 


shape and peak of the SED and the optical depth in the 9.7 pm 
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Wavelength (/xm) 


Fig. 3. The left column shows the density distributions (grey scale) and the location of forsterite (green) for the different models described in 
sect. |4. 1| The white arrow shows the inclination of the system. The second and third columns show the strength of the 11.3 /rm band and the other 
mid-infrared bands and the fourth column shows the flux in the 69 band of forsterite, all as described in sect. |4.3| The different colours of the 
curves in the second, third and fourth column are for the observations (grey) and models with forsterite grains with different grain sizes: black 
(0.1 /rm), blue (1.0 /jm), red (2.0 /rm), orange (4.0 /rm), yellow (6.0 yum), purple (8.0 yum). The forsterite abundances are listed in Table|^ 
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absorption band of amorphous silicates, in addition to reproduc¬ 
ing the resolved images. The best fit of Meixner et al.| ( |2002 1 
required amorphous grains with an exponential cut-off at 200 
pm in size in order to reproduce the mm fluxes, as well as small 
amorphous grains (down to 0.001 pm) to fit the near and mid- 
infrared part of the SED properly. For the central star we use a 
temperature of 5200 K. As shown by [Sanchez Contreras et al.| 
( 2008| l this temperature might be too low. Changes in the tem¬ 
perature of the central object will change the dust temperature in 
the inner region of the circumstellar environment. However, as 
shown in sect. 4.4 and by de Vries et al. ( |2012 1 , this will have no 
significant effects on the modelling of the dust in the case of very 
high optical-depth environments, because then the inner parts of 
the outflow are not directly observed 

We have reproduced the best-fit model of [Meixner et al. 
([2002j l using the radiative transport code MCMax (Min et al. 
[2009 [ l. The MCMax code has been widely and successfully ap¬ 
plied to model observables in a variety of environments (Mul 


[ders et ar]|2011[ [Min et aLl[2013[ [Lombaert et al.[[2012| [201 3| l. 

We refer to these papers for a full description of the features of 
MCMax. MCMax first computes the radiative equilibrium tem¬ 
perature stratification throughout the circumstellar environment. 
The temperature structure is calculated for the different grains 
present in the circumstellar environment separately, so no ther¬ 
mal coupling between different types of dust grains is assumed. 
For the size distribution of the amorphous grains, the opacity 
of grains with different sizes are weighted and combined. Using 
the resulting temperature structure, spectra can be obtained by 
ray-tracing. Dust scattering can be treated in an angle-dependent 
way, or, to speed up the computations ( Min et al.|2009| l by using 
isotropic scattering. In this work we use isotropic scattering. 

We use the same size distribution and optical constants as 
[Meixner et al.[ ( [2002 1 for the amorphous olivine (MgFeSi 04 ), 
namely those of Dorschner et al. (19951, and we extend the 
model of [Meixner et al.| l 2002[ ) by including forsterite grains. 
MCMax uses temperature-dependent opacities, meaning the 
model correctly predicts the 69 pm band shape depending on the 
temperature of the forsterite grains (Mulders et al. 2011| l. For 
crystalline olivine we use the optical constants of forsterite from 
Suto et al. ( [2006 1 , which are available for grain temperatures of 


50, 100, 150, 200 and 295 K. 

The opacities of the dust grains are calculated using the Dis¬ 
tribution of Hollow Spheres ( [Min et aL] ( [2003| l, fmax=10). Other 
particle shape models used in the literature are spherical grains 
and the Continuous Distribution of Ellipsoids (CDF, [Bohren & 
[Huffman! 1983| l. We do not consider spherical grains because it 
has been shown that for crystalline dust species this model pre¬ 
dicts spectral features with shapes and wavelength positions that 
do not compare with spectra of any astronomical object ( [Min[ 
[et al.[[2005| ). [Min et ah] ( [2005 [ l also show that the DHS model 
reproduces the properties of forsterite spectral bands in spectra 
of astronomical objects very well. We also do not consider CDF 
grains for two reasons. The first is that it predicts very similar 
spectral features in shape and position as DHS (see Fig. [^. A 
difference between spectral features of CDF compared to 0.1 pm 
DHS grains is that they are slightly weaker. Second, CDF is only 
valid for particles small compared to the wavelength of the radia¬ 
tion considered and does not allow us to model grains of different 
sizes, which is essential for this work. 


4.3. Extracting the forsterite spectrai features 

A possible 11.3 pm band would be situated on top of the red 
shoulder of the 9.7 pm band of amorphous silicate, which is 
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Fig. 4. The spectral energy distribution of IRAS 17150-3224 in black 
dots. The 8.5 and 10.0 pm points are from jUeta et al.| ( |2000| l, the 12, 
25, 60, and 100 pm points from IRAS, the 800.0 and 1100.0 pm points 
from [Hu et al.[ l [T993[ l and the 1.3 mm from lGuertler et aT|!1996| l. In 
grey is our model based on the parameters of [Meixner et al.[(|2002|). 


in absorption (see the inset of Fig. [T]). The 11.3 pm band has 


been extensively studied in massive AGB stars by de Vries et al. 
( |2010| l and [de Vries et al.[ ( [2014[ l. They study the ratio of the op¬ 
tical depth of forsterite relative to that of the amorphous contin¬ 
uum dust in the 9.7 pm band. This ratio (here called the strength 
S of the 11.3 pm band) is interesting because in the optically 
thick limit it is proportional to the opacities (and abundance) of 
forsterite. 


'Tfo Afo ■ A'fo 
= - OC - 

"^all ^all 


( 2 ) 


where Tof and Taii are the optical depth of forsterite only and the 
optical depth of all dust species together, respectively, and Afo 
is the abundance of forsterite. Since the circumsteller environ¬ 
ment of IRAS 17150-3224 is very dense and the optical depth 
is much larger than one, equation]^ holds. This tells us that if 
there is forsterite in the line of sight to the star and the crys¬ 
tals do not approach a grain size of one micron, an absorption 
band of forsterite must be seen on the red shoulder of the 9.7 pm 
band of amorphous olivine. If the forsterite grains are a micron 
or larger in size the strength S will quickly go to zero. 


We briefly describe how we measure the strength of the 
11.3 pm band (eq. ^ and refer to [de Vries et al.[ ( [2010| l for fur¬ 
ther details. In order to calculate the t^h and we first con¬ 
struct a continuum over the 9.7 pm band of amorphous olivine 
(as shown in Fig.[2l. The optical depth in the 9.7 pm band is cal¬ 
culated using T = -ln(T'/Fcont). where Fcont is the flux level of 
the continuum constracted over the 9.7 pm band of amorphous 
silicate. The optical depth in the 9.7 pm band due to forsterite 
(Tfo) and due to the other dust species (Tou) is now determined 
by constructing a continuum in r-space, under the 11.3 pm band 
(s ee again Fig.[T]l. F or the 69 pm band, we follow the approach 
of Blommaert et al. ( [2014 1 by fitting the continuum with a lin¬ 
ear function (see Fig. [ij. The mid-infrared bands are extracted 
in the region of 21 to 37 pm by fitting a spline to the continuum 
points. The continuum points are taken at (all in pm) 21.7-21.9; 
26.5-26.7; 29.4-31.8; 36.7-37.8. 


10 




















































































B.L. de Vries et al.: Micron-sized forsterite grains in the outflow of IRAS 17150 



Fig. 5. Temperature gradient of the amorphous silicate (solid line) and 
the forsterite (dashed line) grains for the Whole model (see sect. |4. l| l. 
Black lines give the temperature structure over the equator, while the 
grey gives it at the pole. The grey vertical dashed line indicates the 
transition between the low-density AGB wind and the superwind at the 
radius of Rsw = 1944 AU. 


4.4. Base model and testing the results of Meixner et al. 
" 200 ^ 


We find that our model SED, like the model SED of |Meixn^ 
jet al.| ( |2002| l, matches the observed photometry well enough to 
study the spectral bands of forsterite (see Eig. [^. Similar to 
the model of Meixner et al.| ( 2002) 1, we slightly over-predict the 
flux at 100 gm, but due to the large amorphous grains added by 
Meixner et al. ( |2002| l the fluxes around 1.0 mm are well repro¬ 
duced. It is outside the scope of this work to try and improve the 
fit to the SED. 

The spectral features we obtain by adding forsterite to the 
model with the best-fit parameters of [Meixner et al.| ( |2002[ ) are 
shown in Eig. We see that the model reproduces the observed 
69 gm band of IRAS 17150-3224 well in strength, but its wave¬ 
length position is too blue and too narrow. This means that for 
this model the forsterite grains are too cold and/or too small. The 
mid-infrared bands of forsterite in the model spectrum are almost 
completely absent, similar to what is observed. The strength of 
the 11.3 fj.m band on the other hand is too high in our model 
compared to the absence of any 11.3 /urn band in the observed 
spectrum. In sectionj^we will show what happens to the spectral 
features of forsterite if the forsterite grain size and geometrical 
distribution are varied. 

Since different values for the total dust mass in the out¬ 
flow and the central temperature of the star have been reported 
( [Sanchez Contreras et al.|2008| [Guertler et al.|199^ we want to 
test if a difference in these parameters has an effect on the spec¬ 
tral features of forsterite in the spectrum of IRAS 17150-3224. 
We investigate this by varying the central temperature and dust 
mass by +10%. We found this has no significant effect on the 
spectral features. We therefor conclude that the fit of [Meixn^ 
et al. ( 2002| l and the parameters they obtained are of sufficient 
quality to reproduce the optical depth and temperature structure 
of the circumstellar environment in order for us to study the crys¬ 
talline dust component in the outflow. 

In Eig. 1^ we show the equatorial and polar temperature gra¬ 
dient for both the amorphous dust and forsterite that we obtain 
with our model. We find that the temperature of the amorphous 


dust component, at the inner radius of the dust shell, reaches 


~220 K, which is the same value as found by Meixner et al. 
P002[ l. The temperature gradients also show that at the equator 
and within the superwind (Rsw < 1944 AU) the optical depth 
is so high that the temperature of the dust is determined by the 
local radiation field of the dust. This means that both the amor¬ 
phous dust and the forsterite have the same temperature. Out¬ 
side the superwind the optical depth drops and the temperature 
structure of the amorphous dust and the forsterite starts to dif¬ 
fer and the forsterite becomes slightly warmer than the amor¬ 
phous dust. The forsterite is warmer than the amorphous dust 
because 1) the amorphous dust can cool more efficiently than 
forsterite since it has higher far-infrared opacities and 2) since 
both dust species are now heated by the radiation field of the 
inner dust shell. This shell radiates most of its energy in the mid- 
infrared, where forsterite has a higher overall absorption opacity. 
The optical depth at the pole is lower than at the equator and the 
forsterite is slightly cooler than the amorphous dust. 


5. Results 

The mid-infrared and 69 gm bands of forsterite we find for the 
observations and models are shown in Eig. In this section we 
will discuss the forsterite features of the five different models 
introduced in section [4T| 

Whole model 

In Eig.j^it can be seen that the 11.3 gm band in the Whole mod¬ 
els for forsterite grain sizes of <0.1 and 1.0 gm are too strong. 
Only when the grain size is increased to >2 gm does the 11.3 gm 
band become weak enough and compare well with the observed 
absence of the 11.3 gm band. The mid-infrared spectral bands 
do not strongly show in the models for any forsterite grain size. 
The 69 yum band for models with forsterite grain sizes below 6 
fim have a wavelength position that is bluer and a shape that is 
more narrow than the observed band. The model with a forsterite 
grain size of 6.0 fim reproduces the observed 69 gm band well, 
while the 8 gm forsterite grain size model has a too broad 69 gm 
band. 

Superwind model 

Similar to the Whole model, the 11.3 gm bands of the Super¬ 
wind models are only weak enough for the models with forsterite 
grain sizes of >2 pm. The Superwind models with grain sizes of 
<2 pm show absorption bands in the mid-infrared, especially at 
~ 23.5 pm. Only models with forsterite grain sizes of >4 pm 
have weak enough mid-infrared bands to reproduce the observed 
spectrum. The 69 pm bands of the Superwind models reproduce 
the observed band very well, except for the 69 pm band in the 
model with 8 pm sized forsterite grains. This model has a 69 pm 
band that is too broad. 

AGB-wind model 

Eor this model the forsterite is so cold that no spectral features 
are seen at 11.3 pm and in the mid-infrared for any forsterite 
grain size. Because the forsterite is so cold, the 69 pm band is 
also very narrow and blue shifted and the models with forsterite 
grain sizes of <4 pm do not reproduce the 69 pm band well. 
At 8 pm forsterite grain sizes the 69 pm band is too broad, so 
only the 6 pm forsterite grain size model reproduces the 69 pm 
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band well. The jagged feature on the red shoulder of the 69 pm 
band for the 8 pm forsterite grain size model occurs because the 
69 pm band develops a mild double-peaked shape at grain sizes 
of >8 pm and low temperatures (Fig.|^. 

Torus model 

Similar to the Whole and Superwind model the Torus model with 
forsterite grain sizes of <2 pm have a too strong 11.3 pm band. 
In the mid-infrared the Torus models show absorption bands 
(especially at 23.5 pm) for all the forsterite grain sizes, except 
possibly for the 8 pm model where the absorption features are 
weak enough to reproduce the observed absence of mid-infrared 
bands. All but the 6 and 8 pm forsterite grain size models repro¬ 
duce the 69 pm band well for the Torus model. 

Disk model 

We had to increase the forsterite abundance for this model to 
90% in order to see a signal at 69 pm. It can be seen that even at 
this abundance the strength of the 69 pm band cannot be repro¬ 
duced. 



A ifim) 


Fig. 6. Shown in grey is the observed, continuum subtracted 69 pm band 
and over plotted in black a radiative transport model {Whole model, 
forsterite grain size of 4 pm). Shifted to the horizontal grey line are 
shown the scaled absorption opacities of ortho-enstatite ( [Chihara et al.| 
|2002) in black. Ortho-enstatite has a resonance around 68.5 pm, indi¬ 
cated with a dashed black vertical line. 


Blue shoulder on the 69 pm band 


The centre and red shoulder of the 69 pm band of 
IRAS 17150-3224 can be well reproduced by some of our mod¬ 
els, while the small shoulder on the blue side is not (see also 
Fig. 1^. The blue shoulder being located at roughly 68.5 pm 
means it is too far to the blue side for it to be caused by a cold 
forsterite component ( Suto et ar]|2006[ Koike et al. 2003| l. We 
suspect that this blue shoulder might be responsible for what 
caused [Blommaert et al.| ( |2014[ ) to note that the 69 pm band is 
too broad for its central wavelength position. An explanation can 
be the presence of ortho-enstatite (MgSiOa, [Chihara et al.|20d2j l, 
which has a resonance at 68.5 pm. See Fig. ^for a comparison 
of the blue shoulder and the opacities of ortho-enstatite. Ortho- 
enstatite has another resonance at 72.5 pm, but this is in the mid¬ 
dle of the two Herschel/PACS bands, making it difficult to study 
the presence of this spectral feature. Several resonances also oc¬ 
cur in the mid-infrared, but they are not detected in the ISO- 
SWS spectrum. However, if the ortho-enstatite grains are of the 
same size as the forsterite grains, the mid-infrared bands would 
be suppressed as well. Since no other resonances are known to 
us around 68.5 pm, we find it likely that ortho-enstatite is present 
in the outflow of IRAS 17150-3224 with a grain size of a few 
micrometres. 


Summary 

We can exclude the Disk model since we cannot reproduce the 69 
pm band with an abundance as high as 90%. We can also exclude 
forsterite grains of size <2 pm for all models. For the Whole, 
Superwind and the Torus model this is based on the strength of 
the 11.3 pm band. For the AGB-wind model this is based on the 
wavelength position of the 69 pm band. 

We can exclude forsterite grain sizes of >8 pm because for all 
models the 69 pm band becomes too broad. It is difficult to find 
a Torus model that does not have strong mid-infrared bands and 
we therefore find the Torus model unlikely. For the Whole, Su¬ 
perwind and AGB-wind case, satisfactory models can be found 
for forsterite grain sizes between 2 pm and 6 pm. 

Article number, page 8 of 



Fig. 7. Ratio of the radiation pressure over the gravitational force (6) for 
forsterite grains of different grain sizes. 


6. Discussion 


All of our models indicate micron-sized forsterite grains not 
smaller than 2 pm and not larger than 6 pm. This shows that the 
forsterite is still significantly smaller than the largest amorphous 
grains (~200 pm). From dust condensation theory this difference 
can be understood, since the formation and growth of crystalline 
solids requires more specific conditions ([Tielens et ah 1998[ Gail 


|& Sedlmayr|1999t jSogawa & Kozasa|1999| l. A high temperature 

(above ~ 1000 K) and pressure together with sufficient cooling 
times are needed to form a crystalline solid. Otherwise an amor¬ 
phous solid is formed instead. In the outflows of evolved stars 
conditions for crystalline dust formation are only met in some 
cases, close to the central star in the inner part of the outflow. 
Because of this it can be understood that only small amounts 
of small crystals can be formed. This in contrast to amorphous 
solids, which can be formed and grown in larger regions of the 
outflow. 

The fact that de Vries et al. ( 2014| l can model consistently 
the 11.3 pm and 69 pm bands in massive AGB stars, means 
the forsterite crystals in the superwinds of these objects are sub¬ 
micron in size. We then find a large difference in forsterite grain 
size in superwinds of AGB stars and in the outflow of the PPN 
IRAS 17150-3224. At this moment there are no reports of large 
amorphous silicate grains in superwinds and the multitude of 
SED fits to superwind spectra in the far-infrared should have 
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found large grains if they were present Pusttanont et al.||2006| 
2013|l. For the massive AGB star OH 26.5-H0.6, [Groenewegen 


( 2012| l finds an amorphous grain size between 0.15 pm and 0.25 
pm. 

The high mm-fluxes of IRAS 17150-3224 together with our 
analysis of the forsterite spectral features show that there is an in¬ 
crease in grain size for both the amorphous and crystalline com¬ 
ponent in-between the superwind and the PPN phase. If there 
indeed exists a hyperwind, the conditions during it must lead to 
very efficient grain growth and crystallisation. We associate a 
very high density with the hyperwind and indeed dust conden¬ 
sation theory suggests that the density is critical in grain growth 
and the formation of cr ystalline silicates (|Tielens et al.||1998[ 
Gail & Sedlmayr||l999[ Sogawa & Kozasa||1999| l. Besides the 


density, the temperature is also critical for the formation of crys¬ 
talline solids. For olivine the temperature needs to be higher than 
its glass-temperature of ~1000 K, for it to condense into crys¬ 
talline forsterite. Since we see no indication of a disk-like struc¬ 
ture, the formation of the large forsterite crystals we see around 
IRAS 17150-3224 likely happened in the hot dust condensation 
region close to the star during the hyperwind, although we can¬ 
not exclude the AGB-wind model. 

The firm lower-limit we find for the forsterite grain size is a 
challenge to explain. If the formation of crystalline material is 
efficient during the hyperwind, why is the forsterite grain size 
dominated by micron-sized crystals and why is there no popu- 


the area to volume ratio of the grains. This shows that the lower- 
limit we find for the forsterite grains around IRAS 17150-3224 
is probably not directly linked to the wind driving by scatter¬ 
ing on micron-sized grains. One could imagine that forsterite is 
the wind driving dust species and that a population of ~0.5 pm 
grains are blown out and that those subsequently grow to the ob¬ 
served size of 2-6 pm, but since the density quickly drops after 
the onset of the outflow, this would require very special circum¬ 
stances. 

Another possible scenario to explain the lower-limit of the 
forsterite grain size is that an initial population of forsterite 
grains is formed with a size distribution ranging from sub¬ 
micron grains up to ~6 pm, but that the smaller grains are pref¬ 
erentially destroyed or amorphisised. If the smaller forsterite 
grains are preferentially destroyed, one has to explain why the 
small amorphous grains are not destroyed along with them. In¬ 
stead of a destruction mechanism for both the amorphous and 
forsterite dust one could consider the selective amorphisation of 
forsterite grains instead. 

One possible amorphisation mechanism is the destruction of 


ample Ar^^ or Fe^^, 

Borg et al. (19801; 

Day (19771; Bradley 

(1994: 

Demyk et al. 

(|2001|l; ICarrez et a 

1. (|20021l; Jager et al. 

(2003: 

Brucato et al.|( 2004P; Kemper et a 

.(2004 12005 l).Kem- 

per et r 

il.|(|2004|l shows that for the intersh 

;llar medium (ISM) it 


takes several million years to destroy several percentage points 


lation of sub-micron crystals? One possible mechanism we ex¬ 
plored is based on the wind driving mechanism in oxygen-rich 
AGB outflows. The outflow of AGB winds is thought to be ini¬ 
tiated by pulsations, which pushes gas to distances from the star 
where it can condense into dust grains (|Wickramasinghe et al. 

by mass of sub-micron forsterite grains. Since the super- and hy¬ 
perwind of IRAS 17150-3224 combined would not even take 
more than several thousand years, these amorphisation mecha¬ 
nisms are too slow. Possibly, more efficient amorphisation can 
be obtained at the inner-rim of the AGB (hyper) wind, where 
the fast stellar wind will impact the recently stopped hyperwind 

1966 

Gehrz & Woolf|| 1971 [ |Sedlmayr|| 19941 |Habing & Olofs- 

son||20031l. Then radiation pressure on the most refractory dust (|Balick & Frank|2002 1 , but that could not easily explain the ah- 


grains drives the outflow to escape velocities. For carbon-rich 
AGB stars this works well ([Winters et al.|2000[ Gautschy-Loidl| 


|et al.|20()4| Nowotny et al.pOlO 201 1[ Sacuto et al.|2011[ ), but 

driving a wind by absorption of photons in an oxygen-rich en¬ 
vironment is difficult (e.g. |Jeong et al.||2003| ). This is because 
the refractory dust grains in an oxygen-rich environment (e.g. 
olivine and alumina) are transparent in the optical, meaning they 
absorb few photons, and thus momentum, from the radiation 
field of the central star (which has a temperature of ~3000 K). A 
suggestion to drive winds around oxygen-rich AGBs that seems 
to work well is momentum transfer to the grains by scattering 
( [H6fner|[2008| |Bladh & H6fner|[2012] l. For this to be efficient, 
the grains need to grow to a size of ~0.5 pm. Since forsterite 
is a refractory solid, it could be a driver of the outflow of mas¬ 
sive AGB stars. We tested if the grain size dependence of the 
scattering efficiency of forsterite grains might be a way of grain- 
size-filtering, which might explain the lower limit for the grain 
size of forsterite around IRAS 17150-3224. 

The radiative transport code MCMax that we use for the 
models (see sect. |4.2| i lets us calculate the p of grains. This p is 
the ratio of radiation pressure to gravitational force on the grain. 
The code takes into account the anisotropic scattering proper¬ 
ties of the grains. A discussion on the driving of a superwind or 
a hyperwind is outside the scope of this paper, but we simply 
want to see how the p of forsterite behaves as a function of grain 
size. Therefor we calculated the p for forsterite grains at different 
grain sizes, shown in Fig. [7] We see, as shown before by Hofner 


( |2008| l, that the momentum transfer by scattering becomes more 
efficient when the grains approach a size of ~0.5 pm, but a max¬ 
imum p is also reached at 0.5 pm and the scattering becomes 
less efficient at larger grain sizes. This is due to the decrease in 


sence of small forsterite grains throughout the whole super and 
hyper outflow. 


7. Conclusions 

From studying the spectral features of forsterite in the spectrum 
of IRAS 17150-3224 we can draw the following conclusions; 

- The forsterite crystals in the outflow are dominated by a grain 
size population between 2 pm and 6 pm. 

- For the forsterite component we cannot distinguish between 
the Whole, Superwind or AGB-wind models. 

- We exclude the Torus and Disk models. 

- The lower-limit for the forsterite grain size is likely unre¬ 
lated to any wind driving mechanism based on scattering on 
micron-sized dust grains 

- The blue shoulder of the 69 pm band cannot be explained 
by cold forsterite. We suggest the possibility of micron-sized 
ortho-enstatite as a carrier. 

We also speculate that the following processes might be re¬ 
lated to the hyperwind; 

- The amorphous grains grow to several times 100 pm in size 
during the hyperwind 

- Forsterite grains can grow up to sizes as large as ~6 pm dur¬ 
ing the hyperwind 

- The absence of sub-micron forsterite grains could be due to 
an amorphisation mechanism currently unknown 
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